鼓泡流化床中流动特性的多尺度数值模拟

鼓泡流化床因其较高的传热特性以及较好的相间接触已经被广泛应用于工业生产中,而对鼓泡流态化气固流动特性的充分认知是鼓泡流化床设计的关键.在鼓泡流化床中,气泡相和乳化相的同时存在使得床中呈现非均匀流动结构,而这种非均匀结构给鼓泡流化床的数值模拟造成了很大的误差.基于此,以气泡作为介尺度结构,建立了多尺度曳力消耗能量最小的稳定性条件,构建了适用于鼓泡流化床的多尺度气固相间曳力模型.结合双流体模型,对A类和B类颗粒的鼓泡流化床中气固流动特性进行了模拟研究,分析了气泡速度、气泡直径等参数的变化规律.研究表明,与传统的曳力模型相比,考虑气泡影响的多尺度气固相间曳力模型给出的曳力系数与颗粒浓度的关系是一条分布带,建立了控制体内曳力系数与局部结构参数之间的关系.通过模拟得到的颗粒浓度和速度与实验的比较可以发现,考虑气泡影响的多尺度曳力模型可以更好地再现实验结果.通过A类和B类颗粒的鼓泡床模拟研究发现,A类颗粒的鼓泡床模拟受多尺度曳力模型的影响更为显著.      

SPH-based simulation of granular collapse on an inclined bed

The paper presents a numerical model for simulating a granular flow and its deposition on an inclined bed. A granular material is described as an elastic–plastic continuum and its constitutive law, namely Hooke's law, is discretized on the basis of the Smoothed Particle Hydrodynamics (SPH) method. In the equation of motion, however, the artificial viscosity, which is widely used in SPH, is not applied. The diffusive term derived from Hooke's law is introduced with a diffusion coefficient that varies depending on the stress and strain rate based on the Drucker–Prager yield function. The model is verified and validated through two numerical tests. It is shown that the basic elastic–perfectly plastic characteristics are reproduced with a simple shearing test. The effects of the diffusion coefficient and spatial resolution are investigated to show the validity of the model. In the simulation of the gravitational collapse of a granular column on an inclined bed, the performance of the model from the final deposition profile, the time history of the front position of the granular flow, the maximum runout distance, and the velocity profile are investigated for several cases of basal inclinations. The calculated results show good agreement with the experimental results.     

Numerical simulation of counter-current flow field in the downcomer of a liquid–solid circulating fluidized bed

A comprehensive study on the hydrodynamics in the downcomer of a liquid–solid circulating fluidized bed (LSCFB) is crucial in the control and optimization of the extraction process using an ion exchange LSCFB. A computational fluid dynamics model is proposed in this study to simulate the counter-current two-phase flow in the downcomer of the LSCFB. The model is based on the Eulerian–Eulerian approach incorporating the kinetic theory of granular flow. The predicted results agree well with our earlier experimental data. Furthermore, it is shown that the bed expansion of the particles in the downcomer is directly affected by the superficial liquid velocity in downcomer and solids circulation rate. The model also predicts the residence time of solid particles in the downcomer using a pulse technique. It is demonstrated that the increase in the superficial liquid velocity decreases the solids dispersion in the downcomer of the LSCFB.     

Estimation of particle dynamics in 2-D fluidized beds using particle tracking velocimetry

The experimental characterization of particle dynamics in fluidized beds is of great importance in fostering an understanding of solid phase motion and its effect on particle properties in granulation processes. Commonly used techniques such as particle image velocimetry rely on the cross-correlation of illumination intensity and averaging procedures. It is not possible to obtain single particle velocities with such techniques. Moreover, the estimated velocities may not accurately represent the local particle velocities in regions with high velocity gradients. Consequently, there is a need for devices and methods that are capable of acquiring individual particle velocities. This paper describes how particle tracking velocimetry can be adapted to dense particulate flows. The approach presented in this paper couples high-speed imaging with an innovative segmentation algorithm for particle detection, and employs the Voronoi method to solve the assignment problem usually encountered in densely seeded flows. Lagrangian particle tracks are obtained as primary information, and these serve as the basis for calculating sophisticated quantities such as the solid-phase flow field, granular temperature, and solid volume fraction. We show that the consistency of individual trajectories is sufficient to recognize collision events.     

A solids mixing rate correlation for small scale fluidized beds

A new first degree solids mixing rate is proposed to evaluate the mixing of solids in small scale fluidized beds. Particle mixing experiments were carried out in a 2D fluidized bed with a cross-section of 0.02 m × 0.2 m and a height of 1 m. White and black particles with average diameters of 850 and 450 μm were used in our experiments. Image processing was used to measure the concentration of the tracers at different times. The effects of four representative operating parameters (superficial gas velocity, ratio of tracer particles to bed particles, tracer particle position, and particle size) on mixing are discussed with reference to the mixing index. We found that the Lacey index depends on the concentration of the tracers. The position of the tracers affects the initial mixing rate but not the final degree of mixing. However, the new mixing rate equation does not depend on the initial configuration of the particles because this situation is considered to be the initial condition. Using the data obtained in this work and that found in literature, an empirical correlation is proposed to evaluate the mixing rate constant as a function of dimensionless numbers (Archimedes, Reynolds, and Froude) in small scale fluidized beds. This correlation allows for an estimation of the mixing rate under different operating conditions and for the detection of the end point and/or the time of mixing.     

Experimental investigation on the flame front resistance of gas channel growth with melt formation in iron ore sinter beds

The resistance of the flame front within the solid bed constitutes a fundamental and crucial area in porous bed combustion as the flame front propagation is highly related to the productivity and product quality. This paper focuses on the iron ore sintering, a thermal agglomeration process in steel mills. The results from a detailed experimental study of the pilot-scale pot tests under the conditions of a wide range of fuel rate are presented. The primary objective is to provide better understanding of the growth of gas channels relating to melt formation in the flame front and its resistance to flow. The sintering bed was divided into several zones based on the temperature profile and component distribution. Even though there is a continuous one-to-one replacement of humidified zone with porous sintered zone, a constant air flow rate during sintering could be obtained, indicating the ～100?mm high-temperature zone has a controlling effect on sintering bed permeability. The specific pressure drop value in high-temperature zone increases from ～3?kPa in upper bed to ～7?kPa in bottom bed, which varies with the bed temperature and structure properties. Both the green bed and sintered bed were scanned by X-ray computed tomography, the reconstruction and image analysis showed that the sintered bed has large gas channels and many more closed pores due to solid-melt-gas coalescence. More melt is generated when the heat is accumulated along the bed or input higher coke content, showing a propensity to suppress the gas channel growth and amplify the mismatch of gas transportation along the bed. Higher coke rate leads to a higher resistance in flame front, resulting in a slower flame front speed. These results are aimed to provide quantitative validation for improvements of a numerical sintering model in a future work.     

Particle mixing rates using the two-fluid model

In this work, a new methodology is introduced to calculate the solids mixing rate in dense gas-fluidized beds using the two-fluid model. The implementation of this methodology into an existing two-fluid model code was carefully verified. The solids phase continuity equation was satisfied using our method, and the sensitivity of the computational results to the time step, computational cell size, and discretization scheme was investigated to determine the optimal simulation settings. Using these simulation settings, the degree of solids mixing was observed to rapidly (exponentially) increase with increasing operating pressure and linearly decrease with increasing bed diameter. Our novel methodology can be applied to analyze mixing processes in large lab-scale beds as an alternative to existing time-consuming simulation techniques such as computational fluid dynamics combined with the discrete element model.     

氟添加方式对NiMo催化剂在浆态床上加氢脱硫性能的影响

以F为助剂,采用完全液相法制备了氟改性的NiMo/TiO₂-Al₂O₃浆状催化剂,考察了氟的添加方式对催化剂在浆态床上4, 6-二甲基二苯并噻吩(4, 6-DMDBT)加氢脱硫性能的影响.采用X射线衍射(XRD)谱、氢气程序升温还原(H₂-TPR)、N₂吸脱附实验(BET)、X射线光电子能谱(XPS)、高分辨率透射电镜(HRTEM)对催化剂进行了表征.结果表明,在不添加浓硝酸的情况下,在活性组分加入之前引入氟能显著提高催化剂的比表面积和孔径,促进反应过程中金属镍在催化剂表面的分散,更大程度地减弱催化剂中活性金属与载体间的相互作用,有效提高Mo的硫化度及MoS₂的堆积层数,从而生成了较多的II类Ni-Mo-S活性相,促进芳香环的加氢和C-S键的氢解,对4, 6-DMDBT有较高的加氢脱硫活性.     

CFD modeling using heterogeneous reaction kinetics for catalytic dehydrogenation syngas reactions in a fixed-bed reactor

A comprehensive 2D computational fluid dynamics （CFD） model was developed to simulate the flow behavior and catalytic dehydrogenation reaction of syngas in a heterogenous fixed-bed reactor （FBR）. The model combined the porous medium CFD model with a reaction kinetics model. To acquire an accu- rate reaction kinetics model, a comprehensive reaction mechanism was studied for the heterogeneous catalytic dehydrogenation reaction ofsyngas over a supported metal catalyst. Based on the reaction mech- anism and a statistical test, a reliable kinetics model was proposed. The CFD model combined with the above kinetics model was validated with one set of experimental data. The CFD model was also used to predict key reaction variable distributions such as the temperature and the component concentrations in the reactor.     

Hydrodynamic behavior of an internally circulating fluidized bed with tubular gas distributors

To better understand the hydrodynamic behavior of an internally circulating fluidized bed, solids holdup in the down-comer （Eso）, solids circulation rate （Gs） and gas bypassing fraction （from down-comer to riser y~R, and from riser to down-comer yRD） were experimentally studied. The effects of gas velocities in the riser and in the down-comer （UR and UD）, orifice diameter in the draft tube （dor）, and draft tube height （HR） were investigated. Experimental results showed that increase of gas velocities led to increase in Gs and yDR, and slight decrease in yeD. Larger orifice diameter on the draft tube led to higher 8sD, Gs and yDR, but had insignificant influence on YRD. with increasing draft tube height, both Gs and YDR first increased and then decreased, while yRD first decreased and then increased. Proposed correlations for predicting the hydrodynamic parameters agreed reasonably well with experimental values.